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A new result for the irN £ term from a George Washington University/ TRIUMF group analysis of irN 
data is presented. The value £=79±7 MeV was obtained, compared to the canonical value 64±8 MeV 
found by Koch. The difference is explained simply by the PSI pionic hydrogen value for a^-p, the latest 
results for the ttNN coupling constant, and a narrower A resonance. Many systematic effects have 
been investigated, including Coulomb corrections, and database changes, and our results are found 
to be robust. In the standard interpretation, our value of £ implies a nucleon strangeness fraction 
y/2~0.23. The implausibility of such a large strange component suggests that the relationship between 
£ and nucleon strangeness ought to be re-examined. 
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^sO . 1 Introduction 

The ttN sigma term £ has long been a thorn in the side of low energy quantum chromo dynamics 
(QCD) [Q,|2|. The canonical result S=64±8 MeV was obtained by Koch ||||] based on an analysis 
of pre-1980 Trp and tttt scattering data, KH80 @,f|. Gasser, et al. || later developed an alternative 
method of extracting £ which agreed perfectly with Koch when using the same KH80 solution. In 
■ the usual picture, the nucleon strangeness parameter is 
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y/2 = < N ^ > (i) 

D . < N\uu + dd\N > 



The canonical E result yields y = 0.11 ± 0.07, whereby the strange quarks would contribute ~110 
MeV to the nucleon mass, an amount considered too large to be physical in light of results from 
e.g. neutrino scattering [jTj. This "sigma term puzzle" spawned a whole generation of ttN scattering 
experiments that have greatly increased the size and the quality of the scattering database. 

A long-standing prejudice has been that new and better ttN scattering data and an updated 
analysis ultimately would result in a smaller value for y. With the new generation of experiments 
almost all completed, our George Washington University /TRIUMF group has sought to extract 
the £ term as part of our ongoing ttN partial-wave and dispersion relation analysis program, which 
employs the most up-to-date 7riV scattering data in our SAID database ||] . Our main conclusion is 
that contrary to wishful expectation, a thorough analysis of the new data has yielded a larger value, 
£ = 79 ±7 MeV, which can be understood simply in light of the new experimental information. The 
sigma term and our analysis will be summarized briefly. Details can be found e.g. in Refs. f||,||,^], 10|. 



2 The Pion-Nucleon Sigma Term 

The sigma term & measures the nucleon mass shift away from the chiral {m u = rrid = 0) limit, 
thereby parameterizing the explicit breaking of chiral symmetry in QCD due to the non-zero up 
and down quark masses. Models of nucleon structure are required to determine a. The canonical 
result a = 35 ± 5 MeV is due to Gasser |ll| based on SU(2) chiral perturbation theory plus meson 
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Figure 1: Determination of the wNN cou- 
pling constant from the Hiiper dispersion 
relation. The y-intercept gives the coupling 
g 2 /M, and the left (right)-hand side of the 
figure is dominated by n~p (tt + p) data. 
This technique is well suited to determine 
the coupling constant since most system- 
atic effects (e.g. Coulomb corrections) affect 
each side asymmetrically, "pivoting" the 
curves about the intercept, hence greatly 
reducing their effect on g 2 . 



loop corrections. One obtains the strangeness y from 



a(0) = ^— (2) 
1-2/ 



where the theorem of Brown, Pardee, and Peccei [12| relates <r(0) to the isoscalar invariant 7riV 
scattering amplitude D + (v,t) at the "Cheng-Dashen" point v = 0, t = 2m^: 

S = F^D + (0,2ml) 

= a(2ml)+A R (3) 

where 

a(2ml ) = (7(0) + A a (4) 

and ^=92.4 MeV is the pion decay constant, v is the crossing energy variable, and t is the four- 
momentum transfer. The "remainder term" A# is small (<2 MeV J14{|). The nucleon scalar form 
factor a(t) shifts by an amount A<j=15 MeV from t = to t = 2m^, calculated from a tttt dispersion 
relation analysis || and recently confirmed by a chiral perturbation theory calculation [15|. The 
bar over D + indicates that the pseudo-vector Born term has been subtracted. 

The Cheng-Dashen point lies outside the physical ttN scattering region, so the experimental 
D + amplitude must be extrapolated to obtain E. The most reliable extrapolations are based on 
dispersion relation (DR) analyses of the scattering amplitudes B. The Koch result E = 64 ± 8 
MeV was based on hyperbolic dispersion relation calculations ||. More recently, Gasser, et al, 
(GLLS) H developed another dispersion theoretic approach based on forward subtracted 7riV dis- 
persion relations. Expanding D + (t) as a power series in t, the experimental sigma term E can be 
expressed as 

E = F^(4 + 2ml4 1 + ...) (5) 
= F2{d& + 2m*d&) + A D (6) 
= Ed + Ad (7) 

The GLLS, or truncated, sigma term E^ is obtained via the subthreshold coefficients d^o ana - ^oi> 
calculated from the forward subtracted D + and "derivative" E + dispersion relations, respectively. 
They can also be determined from the subtraction constants D + (0,t) = C + (0,i) in the fixed-t 
dispersion relation C + (v, t). The intercept of the curve D + (0,t) yields d^, whereas the slope at 
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t = yields d^. The "curvature correction" term = 12 ± 1 MeV was determined from a 7T7t 
dispersion relation analysis Q. The great advantage of this approach is that <r(0) can be obtained 
simply from via || 

cr(0) = S d - (3 ± 3)MeV (8) 

since the correction terms A CT and almost cancel, both having similar ttit amplitude input ||. 

The analysis of Ref. || used the Karlsruhe KH80 || 7riV phases as input and fit just the low 
energy data. Their result was ~50 MeV, or £ ~62 MeV (with Ao=12 MeV), in agreement with 
Koch p]. Questions regarding the accuracy of the P + dispersion relation integral, which is more 
sensitive to the smaller and more poorly known higher partial waves than other dispersion relations, 
were answered by the good agreement which demonstrated the reliability of the approach. 

3 Analysis Procedure 

Solutions from our ongoing ttN partial-wave and dispersion relation analysis are released when 
changes to the database and analysis method warrant S. Details of our analysis method can be 
found in Ref. JS , 10|, 16 ] . An energy-dependent ttN partial- wave analysis (PWA) is performed on 



the available data up to 2.1 GeV pion laboratory kinetic energy, applying constraints from forward 
C^w) and E ± (uj) DRs, as well as fixed-t B±(v,t) (in the "Hiiper" form f§) and C^i/.t) DRs. 
These dispersion relations are constrained [] to be satisfied to within <2% from 30 to 800 MeV for 
—0.4 < t < 0.0 GeV 2 /c 2 . The dispersion integrals use the Karlsruhe KH80 phases from 2.1 to 4.5 
GeV and high energy parameterizations above that using forms found in Refs. (j, 17|. 



Dispersion relations depend on a priori unknown constants e.g. scattering lengths and g 2 . Our 
analysis determines these constants by a best fit to the data and the dispersion relations. The 
coupling g 2 , the ir~p s-wave scattering length a n - p , and the p-wave scattering volume a± + were 
fixed for each fit over a grid of values (for reasons of fit stability) , where the combination with the 
lowest x 2 yields the final solution. The fitting procedure automatically chooses the best-fit isovector 
scattering length ag + and volume oj, , and the subtraction constants C^(0, t). This method enables 
us to check their sensitivity to various systematic effects, e.g. database changes. 

The low energy Pi 3 partial wave is constrained to follow the expected partial wave dispersion re- 
lation behaviour in its Chew-Low approximation form [l£| , which our other p- waves satisfy without 
constraint. As well, the low energy F and higher partial waves, too small to be determined from the 
irN scattering data, are constrained to agree with those calculated by Koch Jl9| from partial wave 
projections of fixed-t dispersion relations, which are dominated by t-channel (tttt) contributions. 
This ensures that our higher partial waves satisfy analyticity and unitarity requirements. 

4 Results 

Our main results are summarized in Figs, [l] and || and Table |l[ We find for the irNN coupling 
constant g 2 /4-7T = 13.69 ± 0.07 (or f 2 = 0.0757 ± 0.0004), stable in our solutions for many years 
(see [16]). Our coupling constant agrees with most other recent results, in particular the compre- 



hensive AW and 7riV analyses of the Nijmegen group (see Ref. [20 and references cited therein). 



Note that this result is perfectly consistent with both the Goldberger-Treiman discrepancy pl|,|22[ 



and the Dashen-Weinstein sum rule [23], removing a long-standing inconsistency when using the 
older Karlsruhe value |J 14.3±0.3. 

For the s-wave scattering lengths we obtain Sa„- p = 0.261 m" 1 and 3oq + = 0.260 m" 1 , with 

T")S1 

1-2% uncertainties. The ir~p scattering length agrees with the PSI pionic hydrogen result 3a^_ p = 

1 This range has been increased from previous analyses, and in practice the dispersion relations are well satisfied 
somewhat beyond that range due to the energy dependent partial wave forms 
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Figure 2: The amplitude C (0, t) (points) eval- 
uated from fixed-t C + {u,t) dispersion relations. 
A fit (dashed line) yields the subthreshold coef- 
ficients in the table. The solid diagonal line is 
inferred from forward C'+ and E + dispersion re- 
lations, and agrees perfectly with d 00 and d m in 
the table. The curvature terms ( doi,i > 2) im- 
ply Ac >11 MeV, consistent with the canoni- 
cal result 12±1 MeV from Ref. @. The ampli- 
tude is very small as expected at t = m% ( "Adler 
point" ) . The overall consistency tends to support 
our result for the sigma term, £ ~ 79 ± 7 MeV. 



0.2649±0.0024 m, , while the isovector scattering length satisfies the Goldberger-Miyazawa-Oehme 



(GMO) sum rule [21 when using our coupling constant and integral Jgmo = —1-08 ± 0.03 nib -1 . 
The p-wave scattering volume a[ + = 0.133 m~ 3 is consistent with recent analyses of low energy 
data |2f| , as expected since the resonanct P33 partial wave dominates the low energy data and aj, . 

The dispersion relations are very well satisfied up to about 1 GeV, in general much better 
than KH80. From the forward C + (oj) and E + (uj) dispersion relations, we obtained the coefficients 
doo = —1.30 m" 1 and cZqi = 1-19 m~ 3 , in perfect agreement with the results from the slope and 
intercept of the C + (0, t) subtraction constants at t = 0, shown in Fig. The equivalent J01 
result from the E + dispersion relation, with its ~ / 3 sensitivity to partial waves, and the C + (v,t) 
dispersion relation, with its ~ I sensitivity, supports the reliability of our higher partial waves. 

Figure || shows a polynomial fit to C + (0,t) near t = 0, from which Jo2> ^03, and J04 were 
estimated. The c£o2 coefficient is in perfect agreement with the Karlsruhe result [Q], while the sum 
of the higher order terms yield a curvature correction Ajj> > 11 MeV, in agreement with the 7T7t 
dispersion relation result || 12 ± 1 MeV. Moreover, the curve extrapolates to about -4 MeV at the 



"Adler point" (i=m^), consistent with expected corrections to the Adler Consistency Condition [26], 
where it would be identically in the chiral limit. Compatibility with t-channel dispersion relations 
and chiral constraints gives us confidence in the reliability of our subthreshold coefficient results. 

From the above subthreshold coefficients, and the curvature correction from Ref. [||, our result 
for the sigma term is £=79±7 MeV, compared to the Koch value 64±8 MeV ||. Though surprising, 
the result is readily explained by the new experimental information. Table || shows the breakdown 
of T,^ into its dispersion relation terms for both the KH80 solution and our own. With respect to 
KH80 result, the new PSI pionic hydrogen and deuterium scattering length |^4j a^, ~ 0.000 m" 1 , 
which we reproduce, causes a 7 MeV increase. With a irNN coupling constant g 2 /Air ~ 13.7 |l§,[20| 



£<2 increases by 6 MeV 0. It is well known that the KH80 solution overshoots the data on the left 
wing of the A resonance. Our solution fits the available data much better than KH80, resulting in 
a narrower A width. This leads directly to the 3 MeV increases in each of the dispersion integrals 
shown in Table |l]. Consequently there is sound experimental evidence to support our new S result. 

4.1 Systematic Checks 

Perhaps the most important systematic check is the sensitivity of our results to the scattering 
database. Around the A resonance, there is a well known disagreement between the TRIUMF tt^ 



differential cross section [27] and PSI tt total cross section data on the one hand, and the 



2 The above increases were also noted in Ref. |Q 
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Solution 


S d [MeV] = 


"OqI const. 


Born 




"ciji " const. 


Born 




KH80 


50 = 


-7 


+9 


-91 


+352 


-142 


-72 


FA01 


67 = 





+9 


-88 


+351 


-136 


-69 


difference 


17 = 


+ 7 





+3 


-1 


+6 


+3 



Table 1: Comparison of from the Karlsruhe solution KH80 || and our recent solution FA01. The change 
in the C + subtraction constant (aj + ) term, the E + Born term, and both integral terms are consistent with 
expectations from, respectively, pionic atom data the coupling constant g 2 /4ir ~ 13.7 (see Ref. p0|), 

and a narrower A resonance width. Values are rounded. See text for details. 



older CERN results |^9| on the other. Only a small increase 0.07 and 4 MeV was observed in </ 2 /47r 
and respectively, for the "CERN-only" database. In practice, both sets are included in the final 
fitQ. We found that weeding out large \ 2 data sets had little effect on the result. Also, since the 



low energy data are consistent with the PSI pionic atom results [24|, we conclude that there are no 
large systematic effects from reasonable changes to the current scattering database. 

The hadronic amplitudes are corrected for Coulombic effects following the Nordita prescrip- 
tion p0|, supplemented in this analysis at high energies by extended-source Coulomb barrier fac- 
tors The current approach improved the agreement with the PSI pionic atom results over our 
previous Nordita+point-source barrier results; however, neither the coupling constant nor varied 
outside the errors when using point- or extended-source barrier factors exclusively, or the Nordita 
corrections supplemented by either. Moreover, the isospin-violating A resonance is "split" defining 
"hadronic" = "A ++ ", consistent with the Nordita definition, but find no difference to our previous 
approach with "hadronic" = "(A + A ++ )/2", or with no splitting at all. We conclude that there 
are also no large systematic uncertainties from our Coulomb correction scheme. 

The implementation of our dispersion relation constraints was also checked. We found that 
every reasonable form for the high energy amplitudes (>2 GeV) yields virtually identical results. 
Agreement between the forward subtracted and fixed-t unsubtracted dispersion relations is good 
for reasonable constraints (i.e. typical experimental error ~2%), but suffers if they become too 
tight, <0.5%. Constraining the low energy P13 partial wave to follow the Chew-Low form lowered 

by 6 MeV, but once corrected, reasonable deviations caused changes much smaller than our 
error bar. We also had solutions where the low energy F and higher partial waves were not rigor- 
ously constrained to the Koch values IS], and no significant difference was found. Furthermore, 



Olsson [32 1, from a new dispersion relation sum rule, and Kaufmann and Hite [33|, from an interior 
dispersion relation analysis, obtained values for X consistent with our own using an earlier SAID 
solution. Consequently, we are confident in the reliability of our dispersion relation analysis. 

5 Summary 

In summary, we have performed a comprehensive partial wave and dispersion relation analysis 
of the available irN scattering data up to 2.1 GeV that includes several improvements upon prior 
analyses |9|,|i"0fl. For the pion nucleon coupling constant we obtained <? 2 /47r = 13.69+0.07, consistent 
with our previous determinations [||,[n| and the Nijmegen results [ p0|. Our s-wave scattering lengths 
agree with the latest PSI pionic hydrogen and deuterium results |24| . Our ttN sigma term result 
is = 67 + 6 MeV, or X = 79 ± 7 MeV, compared to the canonical result 64 ± 8MeV from 
Koch H . These results have proven robust with respect to the many systematic checks that we have 
performed. In light of the large nucleon strangeness content y/2 ~ 0.23 inferred in the standard 

3 The other low and A resonance energy data are fit somewhat better in "TRIUMF+PSI-only" solution 
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picture, we believe that alternate interpretations of a large sigma term ought to be examined 
carefully. 
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